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22 Mr. Hedrick , Probable Error liii. i, 

similar to his, with a provision for carrying the observer similar 
in general character to his, making it, in fact, exactly on the 
lines he has laid down. 

The first modification would, perhaps, be the best form 
altogether, as the sides and ends of the outer house could be 
made with a considerable slope, so as to reduce the upsetting 
power of the wind. This would retain the important position of 
the observer as close as possible to the centre of motion of the 
telescope. 


Probable Error of the Cloclc Correction when both the Clock Bate 
and the Instrumental Constants are found tyy a Least Squares 
Solution of a Single Night’s Observations. By the Rev. 
John T. Hedrick, S.J. * 

(Communicated by the Bev. W. Sidgreaves, S.J.) 

It seems to be a common practice in the comparison of time 
determinations, as in longitude work, to give a constant probable 
error for the clock correction of each night and to combine the 
different nights* work according to these probable errors. This 
is equivalent to taking the clock rate as known perfectly, whereas 
its value has also its degree of uncertainty, making the probable 
error of the clock correction quite different for different moments. 
Of course, if the time comparison takes place during the period 
of observation, and especially if about its middle, there will be 
less objection to the method. 

A reason for this neglect of the probable error of the clock 
rate may be the difficulty of determining it when another quite 
common practice is followed of determining the clock rate from 
the comparison of the clock corrections found from night to 
night. If results are expected of sufficient accuracy to make it 
worth while taking account of probable errors, this practice 
seems hardly defensible for observations made with portable 
instruments. For we are likely to have periodic or irregular 
fluctuations of rate in a clock, not of the best, exposed te 
receiving injurious shocks in carriage, and not set up so as to be 
thoroughly protected against the effects of change of tempera¬ 
ture, or even against instability in the clock pier. Whatever 
periodic or accidental errors may exist in the star places of the 
various modern ephemerides, they are so small that, except in a 
fixed observatory, the proper method is to find the clock rate 
each night from that night’s observations. 

Supposing that a least squares solution of the night’s work 
is thought worth the while, there is still a difficulty in getting 
the general expression for the probable error of the clock correction 
at any moment. Its probable error at a given moment may be 
found by making this moment the epoch of the constant part of 
the total clock correction; but this may bring in practical diffi.- 
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cnlties in computing, requiring, for example, the use of more 
decimal places in the coefficients and logarithms. 

Let 

AT be the clock correction at any moment, 

AT 0 the clock correction at some assumed moment, and 
p the clock rate. 

From the sum of the squares of the residuals after substitu¬ 
tion in the equations of condition of the values found by the 
least squares solution and the coefficients occurring in the solu¬ 
tion, we can find the probable error of AT 0 and p. The probable 
error of AT will be a function of these probable errors and the 
interval between the epochs of AT and AT 0 , but it will not be a 
function of these only, if the proper epoch is not chosen for 
AT 0 ; for in this case the probable error of AT 0 will itself 
depend on that of the rate. (See equations (2) and (3), below.) 

Let 

n\ .... 

be the known terms of the equations of condition after reduc¬ 
tion to weight 1 . AT 0 and p are each determined as a linear 
function of the ns. 

Let 

AT 0 =aV + aV , + a"V ,, + . . . . 
p = j8W + j8'V' + j 3 ''V'' + .... 

where a', / 3 r , <fcc., are functions of the constant coefficients of the 
equations of condition. 

Now 

AT = aT 0 + t. p, 

where t is the interval between the epochs of AT and AT 0 . 
Hence 

AT=(a' + / 3 ' . *>'+ («" + /8" . *)»" + («'" + j8"' . t)n’" + .... 

Let p be the probable error of an observation of weight 1, 
and^AT, pAT 0 , andjppthe probable errors of these quantities. 
We have (Chauvenet, ii», pp. 498,499, and 516)*— 

(pAT 0 y=p 2 . [cca], {pp) 2 =f . [/ 3 i 8 ], 

(pAT y=f . [(a + j8 . ffta + fi . 0 ] =P 2 . ([«a] + 2[afl . t+ [fi 0 ] . <*). 

From the degree in t of the parenthesis it is evident that there 
is a value of t which makes p AT a minimum, and but one value. 
Let us put AT 00 = AT for this epoch. It is evident also that if 
the epoch of AT 00 were chosen as the epoch of AT 0 , we should 
have [a P] = o; for if it is not zero pAT will not be a 
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24 Mr. Hedrick, Probable Error liil i, 

minimum for that epoch. Hence, counting t from the epoch 
of A T 0 o— 

M + [^].O = (?at 00 )H(pp) 2 .^ ... (I) 

Hence the probable error of AT 00 is independent of that of 
the rate. If the epoch assumed for AT 0 be not that of A T 00 , and 
if be the interval between the two epochs, 

(pATJ^tpAT^ + ^p)*.^.(2) 

and 

(i»AT)*=(|iAT 0 )*-(pp)*.<'« + (pp)*.<*.(3) 

There is no difficulty in assigning the epoch of AT 0 o when 
only it and p enter our equations. It is the mean (simple or 
weighted) of the times of observation. But we have commonly 
to determine also the azimuth (a), and in the case of a portable 
instrument the collimation (c) as well, and perhaps the azimuth 
after reversal ( a '). In this case the epoch of AT 00 is not the 
mean of the times of observation, and it is the object of this note 
to show how to find this epoch from the coefficients formed in solving 
the normal equations for the epoch of AT 0 , or to find the quantity t’ 
of equation (2). The probable error of the clock correction at 
any moment will then be found by (1) or (3). 

Since the squares of the probable errors of the quantities 
found by a least squares solution are inversely proportional to 
the coefficients each would have in the final equation if it were 
the first unknown found, it follows that the moment at which 
the clock correction has the minimum probable error is the 
moment which, chosen as the epoch from which to compute the 
ooefficients of the rate in the equations of condition, would bring 
out AT 0 in the final equation with the maximum coefficient. 

Let us suppose the most general case, namely, that we have 
to determine a, a', c, p, and AT 0 , and let w = weight. We have 
this group of five normal equations, using Gauss’s notation. 

[waa]x + [wab]y + \wac\z + [wad]u + [wae]v + [wan] = o, 

~\w&b]x + [wbb]y + [wbc]z + [wbd]u + [wbe]v + [ wbn ] = o, &c. 


„ Let the order of the unknowns and that of their elimination 
ho that given above. 

For [ee . 4], the coefficient of v(AT 0 ) in the final equation, 
or the quantity to be made a maximum, we have 




In the second member all the terms except the last are inde¬ 
pendent of d, the coefficient of the rate, and hence of the epoch 
of AT 0 , and all are positive. Hence [ee . 4] will be a maximum 
when the epoch is so chosen as to make the last term zero, or, 
since [dd . 3] is not infinite, when [de . 3] = o. 
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We have the following expressions, in determinant nota¬ 
tion *:— 

[de . 3] = §waa][wbb][wcc][wde]\ ~ \[waa][wbb][wcc]\, 

[ ee . 3] = \[wad] [wbb] [wcc] [wee ]| -r | [waa] [wbb] [wcc ]|. 

Or, since e, the coefficient of AT 0 in the equations of con¬ 
dition, = 1— 

[de . 3] = \[waa][wbb\]wcc][wd]\ •*-1 [waa][wbb][wcc§, 

[ee . 3] = | [waa] [wbb] [wcc] [«/]| -f-1 [waa] [wM] [wcc ]|. 


The only coefficient affected by a change of epoch for AT 0 is d, 
the coefficient of the rate. The denominator of \_de . 3] is inde¬ 
pendent of d, and, as it is not infinite, to make [de . 3] = o its 
numerator must become zero. This numerator expressed in the 
fuller determinant form is— 


[waa] [wab] [wac] [wad] 


No. of [de . 3] = 


[wab] [wbb] [wbc] [wbd] 

[wac] [wbc] [wcc] [wed] 


[wa] [wb] [wc] [wd] 


Let us suppose that we have formed and solved our normal 
equations «with 4 he epoch T for AT 0 , and let us change the epoch 
to T' = T-H'. The coefficient of the rate is the time of observa¬ 
tion minus the epoch of AT 0 . We shall have the following rela¬ 
tions between the two sets of coefficients, d representing in both 
sets the coefficient of the rate for the epoch T. 

Epoch T, Epoch T'=T + t' 
d, d—t', 

[wad], [wad] — [wa] . V, 

[wbd], [wbd] — [wb] . t', 

[wed], [wed] — [wc] . V, 

[wd], [wd] — [w] . f, 


*, These expressions, as well as those for the other coefficients occurring in 
the solution of the normal equations (Chauvenet, ii. pp. 531 and 532), can he 
found by expressing each as a determinant of the second order; e.g. 

! [dd . 2 ][de . 2]| 

[ee . 3 ]=*^^—. x 


[dd . 2 ] 


[de . 2] [ee . 2]\ 


substituting in each the determinant expressions for the coefficients of lower 
order, and applying the theorem 


|a n a 2 i 


a n a 3J 


a n a 4i 

| • 

«11 ®W1 

1 



a is a as 


a i4 a u 

1 .... 


1 


-.a™- 2 

11 


1 ^ 11 ^ 22^33 • • • • a v 


The first member is the determinant whose element in the /rth row and 
Ith column is 


a n a (k+i)i 
^UZ+l) <Z(fc+l)(!+l) 


Ki-1,2,3 — (»- 1 ). 
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Let us substitute tbe values for the epoch T' in the numerator 
of \_de .3]. We have 

\waa][wab]\wac]\wadi\ — \wa\ . t’ 
[wab][wbb][wbc][tvbd] -*[wb] . f 
[wac'] [wbc ] [wco] \wcd'\ — [ wc ] . f 
[1 wa ] [• wb ] \wc\ [1 wd ] - [ w ] . t' 

Splitting up this determinant into two and factoring out t r 
in the second, we have, for the epoch T 7 — 

No. of [<fe . 3] = | [waa] \wbl~\ \wcc\ [wd] ]—| \_waa\ \_wbU\ [wcc] [w] |. if. 

Putting this equal to zero and noticing that the first term 
and the factor of t' in the second term are proportional to \de . 3] 
for the epoch T and [ee . 3] respectively (4), we find for t' the 
interval from the epoch assumed for AT 0 to that which would give 
the maximum coefficient for AT 0 in the final equation , or to the 
epoch at which the clock correction has the minimum probable error , 

\de . 3] for the assumed epoch 
~ [ee . 3 ] . 

The value found for t' is the quotient of the same coefficients, 
after the equations have been freed from all the unknowns except 
AT 0 and p, of which it would have been the quotient, had there 
been only these two unknowns in the original equations, namely, of 
the coefficient of AT 0 in the normal equation in p by that in its 
own normal equation. It will, in general, differ from \wde~\ 
[wee] = the mean of the times of observation. 

In the example given below the square brackets are used 
also to indicate the number whose logarithm they enclose. 

Summary, 

The epoch of the clock correction of maximum weight, or 
minimum probable error, is not, in general, the mean of the 
times of observation when, besides the constant clock correction 
and the clock rate, instrumental constants are also determined 
from the observations, (5). 

If these quantities are found by a least squares solution, 
this epoch is before or after the epoch assumed for the constant 
clock correction by an interval which is the quotient of the co¬ 
efficient of the constant clock correction in the normal equation 
for the rate by its coefficient in its own normal equation, after 
the elimination of the other unknowns, (5). 

If we count from this epoch, the probable error of the clock 
correction at any other time is what it would be if the constant 
correction and the rate were independently observed quantities— 


No. of [de , 3] for epoch T' — 
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Nov. 1892. Mr. Christie , On Nova Aurigce. 27 

that is, its square is the sum of the square of the probable error 
at this epoch and the product of the square of the probable 
error of the rate into the square of the interval from this 
epoch, (1). Hence the square of the probable error of the clock 
correction at this epoch is equal to the square of the probable 
error of the clock correction at the assumed epoch minus the 
product of the square of the probable error of the rate into the 
square of the interval .between the two epochs. 

Example. 

Number of stars observed, 20. Number of quantities determined, 5. 

Putting v here for the residuals, \wvv\ = 0*^025 r. 

Por an observation of weight 1 ,jf - o 8 *000761 = [6*8815]. 

Mean of the times of observation = 8^98. Assumed epoch of AT 0 = 9 h -o. 

[ee . 4 ]=3-28, (^AT 0 ) 2 = [6*8815] — [0*5159] = [6*3656] =o 8 *ooo232. 

[dd . 4] =444, (pp) 2 = [6*8815] -r [0-6474] = [6-2341] = o 8 -oooi7i. 

[de . 3]= +575 = [ 0 * 75971 , [ee . 3]= +6-85 = [0-8357], [9*9240] = o b - 84 . 

Epoch of maximum coefficient=9^84. 

(pAT 00 ) 2 = o 8 *000232 — o 8, oooi7i . (o*84) 2 =o b cooiii, t?aT 00 = ±o 8, oii. 

At time T, pAT = ± Vo 8 *oooiii +o 8 -oooi7i . (T-9 h *84) 2 . 

Georgetown College Observatory, 

Georgetown , D.C., TJ.S.A. 


On the Photographic Magnitude of Nova Aurigce , as Determined at 
the Royal Observatory , Greenwich. II. By W. H. M. Christie, 
M.A., F.R.S., Astronomer Royal. 

In the Monthly Notices for last March (vol. lii. p. 357) results 
for the photographic magnitude of Nova Aurigce were given as 
determined up to March 9. Since then photographs have been 
taken up to April 1, when the photographic magnitude had fallen 
below 14, and again at the end of August and beginning of 
September, when the star had brightened. Further measures 
were also made of the earlier photographs, modifying slightly 
the results given in the previous paper. The measures were 
made and reduced in the manner explained in that paper, the 
magnitude of the Nova being inferred by comparison with four 
Argelander stars of 8-9 magnitude by means of the formula 

m = 2*5 (log t — 0*97 Vd) + const. 

given in the Monthly Notices , vol. lii. p. 146, the measures of 
diameter of the four comparison stars being used to determine 
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